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[57] ABSTRACT 

Flow rate measuring method and apparatus for deter- 
mining a flow rate of a fluid flowing through a pipe by 
using a fuzzy inference system which is operative on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 
and membership functions defining meanings of state- 
ments described in said IF and THEN parts. Row ve- 
locity of .the fluid is measured at a plurality of locations 
within the pipe by a plurality of sensors. A table of 
membership functions indicating grades of fitness of the 
flow velocities to a distorted flow model distribution is 
provided in a storage together with a mean flow rate 
calculating expression for deriving a flow rate corre- 
sponding to the distorted flow model distribution. A 
membership grade of the measured flow velocity to the 
distorted flow model distribution is determined by re- 
ferring to the membership function table through the 
medium of the measured flow velocity, and a weighted 
mean of the flow rates given by the mean flow rate 
calculating expression is determined by using as a 
weight the membership grade to the distorted flow 
model distribution. 

24 Claims, 17 Drawing Sheets 
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out the same which can evade the shortcomings and 
METHOD AND SYSTEM FOR MEASURING problems 6f the techniques mentioned above. 
FLUID FLOW RATE BY USING FUZZY Another object of the present invention is to provide 

INFERENCE a method and a system for measuring a flow rate which 

5 are capable of measuring a true mean flow rate even in 
BACKGROUND OF THE INVENTION the presence of distorted flows and disturbances by 

The present invention relates to a method and a sys- 4 inf «ence and a fuzzy integration in deter- 

tem for measuring a flow rate of a fluid flowing within mu »««»> and calculation of a mean flow rate of a fluid 
and through a pipe. More particularly, the invention is flowin8 trough a pipe on the basis of outputs of plural 
concerned with a flow rate measuring method and a 10 sei J sors msta «ed within the pipe, 
system therefor with are capable of determining or In view of the above mA other ^j* 5 " which "® 
measuring stably a mean flow rate of a fluid flowing more apparent as the description proceeds, there is 
through a pipe such as tube, conduit or the like in which P mv } ded according to a first aspect of the present in- 
distorted flows and disturbances such as turbulences * entK>n a A?w rate measurmg method of determining a 
may exist 13 ow rate °* * " uld "Owing through a pipe by using a 

For measurement of flow rate of a fluid flowing fuzzy iMer ? ,CC Sy ? tcm . «»P*w *e *fp of 

within a pipe, there has heretofore been known a plural- ""•«»?« fl ° w of the fluid at a plurality of 

sensor type flow meter system in which a plurahty of ? e pi !*> P"*""* ■ ttble "J 

detection elements (hereinafter referred to as the sen- , ft sh, , p ^ nct ' ons » d, " tm 8 of fitness of the flow 
sor) are disposed within the pipe and a single-sensor 20 ^ £ 1*7 ?* di »^U<>*, pro- 
lyl* flow meter system in which- only one sensor is "° W r " e c ?. col » t,n 8 W"* f ° r den ™8 
installed within the 'pipe. $ "°. w ? te ^Pondmg to the distorted flow model 
As a typical plural-sensor type flow meter, there may ^S^TT" .* T^FfS* ^1 
be mentioned those disclosed in JU-A-61.1954188 and „ K£? L °, 1°™? ? d,st0 ^ d fc flo ^ m °f el 
JP-A-61-120016, respectively. The first mentioned flow 25 ftK??. ft.' T 8 funC ? 0n 

the individual sensors are integrated for determining oJ * _ K „ tVtA „° n ^„, . . 

by using a « culating circuit of a corresponding struc- 30 flow 6 model d fstribution. 

1^ nc°^ n ^ of Presence of distortions or devi- According t0 another upetii of the invcntion| ^ is 

ow v K e °? ty t (l K e - a dlst0 ? ed fow velocity ^ a £ QW rate meas ^ for detcmini 

Sfi?T } f y g !° metncal J*™ < e f • a How rate of a fluid flowing through a pipe by using a 

shapes) of « i upstream pipe section, as is disclosed >n fu22y in f cre nce system which comprises Tplurality of 

t" , . , , 35 flow serisors for measuring flow velocities of the fluid at 

In the smgle-sensor type flow meter, a single hot ware a pIuralit of ]ocations within the { a ^ 

element is installed within the pipe under tension m a first membership function table indicating 

order to determine a curvilinear integral (hue integral) grades of fltness of a mean vaJuc of flow velocities 
of a flow velocity for thereby obtain a mean (average) ^thin the pipe to a plurality of distorted flow velocity 
flow rate, as disclosed m JP-A-54- 125060. 40 levels, respectively, a storage containing a second mem- 

The plural-sensor type flow meter in which a plural- bership funct ion table indicating grades of fitness of the 
ity of sensors are installed within the pipe suffers from a flow velocities to a plurality of distorted flow model 
problem that when any one of the sensors is located in distributions for the flow velocity levels, respectively, a 
a region where flow disturbances due to vortexes, flu- stor age containing mean flow rate calculating expres- 
idal separations or the like brought about by geometry 45 5 i ons for deriving flow rates corresponding to the dis- 
or shape of the pipe are significant, infinitesimal varia- tort ed flow model distributions at the flow velocity 
lions or fluctuations (fluidal noise) of the flow velocity levels, respectively, a unit for determining a mean value 
are directly fetched by the mean flow calculating circuit 0 f the flow velocities measured at the plurality of loca- 
as the input thereto, making it difficult or impossible to tions within the pipe to thereby determine on the basis 
obtain the mean flow rate in a stable manner, because 50 of the mean value first membership grades of the mean 
weights imparted to the outputs of the individual sen- value to the flow levels, respectively, by reference to 
sors are selected to be equal to one another or fixed. the first membership function table through the medium 
The single-sensor type flow meter mentioned above of the mean value, a unit for determining second mem- 
also suffers the similar problem because no measures are bership grades of the flow velocity values to the dis- 
taken against the influence of the infinitesimal variations 55 torted flow model distributions for the flow velocity 
in the flow velocity. . levels, respectively/by reference to the second member- 

It is further noted that the location where the infini- ship function table through the medium of the measured 
tesimal variation takes place in the fluid flow velocity is flow velocity values, and a unit for calculating a 
movable as a function of change in the distorted flow weighted mean of the flow rates given by the mean flow 
distribution profile or pattern. Consequently, the mean 60 rate calculating expressions, respectively, by using as 
flow rate calculation using the fixed weight can not weights the first and second membership grades, for 
stably or consistently assure the availableness of accu- thereby determining the flow rate of the fluid within the 
rate flow rate value by flexibly coping with the infinites- pipe. 

imal variation, giving rise to another problem. The present invention starts from the facts that the 

SUMMARY OF THE INVENTION 65 *™ ****** ^ J". 1 ^ 

* through a pipe exhibit a frequency distribution having a 

It is an object of the present invention to provide a unique mean (average) value and variance due to influ- 

flow rate measuring method and a system for carrying ences of distorted flows and disturbances and that the 
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frequency distribution profile is characteristically inher- 
ent to the geometries of the pipe line and the mean flow 
velocity prevailing within the pipe. 

For calculating the mean flow rate by taking advan- 
tage of the characteristics mentioned above, it is pro- 5 
posed as an embodiment of the invention that the de- 
tected values outputted from a plurality of sensors dis- 
posed within the pipe are made use of as fuzzy variables, 
wherein software and hardware are so structurized that 
correspondences are established by the fuzzy inference 10 
between or among the fuzzy variables, a membership 
function table (first membership grade function table) 
containing information about the levels to which the 
flow velocities belong, a distorted flow model function 
table (second membership grade function table) con- 
taining distorted flow distribution model patterns and 
an empirical expression (mean flow rate calculating 
expression) for determining mean flow rates corre- 
sponding to the flow velocity levels and the distorted 2 q 
flow distribution patterns, respectively. 

When a distorted flow distribution of a given pattern 
exists in the vicinity of plural sensors, the fitness grades 
(similarity grades) or membership grades of the dis- 
torted flow distribution to a plurality of distorted flow 25 
distribution patterns prepared previously independent 
of the former are determined by fuzzy arithmetic opera- 
tion, whereon the mean flow rates which correspond to 
the distorted flow model patterns, respectively, are 
determined in accordance with the empirical expres- 30 
sions corresponding to the distorted flow model pat- 
terns. Further, by using the membership grades (fitness 
grades) as weight, a weighted mean value is calculated 
for the distorted flow model patterns. 

By determining the membership grades for the dis- 35 
torted flow model patterns, as mentioned above, it is 
also possible to calculate mean (average) flow rates for 
those distorted flow distributions which are deviated 
from the previously prepared distorted flow model 
patterns. In this manner, the capability of the flow rate *® 
measuring system to cope with variations of the dis- 
torted flow distribution can significantly be enhanced. 
Further, since the mean flow rate is not straightfor- 
wardly determined by using directly the output values 
of the sensors, the mean flow rate as calculated is essen- 45 
tially insusceptible to direct influence of local disturb- 
ances taking place around particular sensors. Thus, 
there can be realized a flow rate measuring system 
which is highly stabilized against disturbances of the 
fluid flow under measurement. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram showing schematically a 
general arrangement of a flow rate measuring system 55 
according to a first embodiment of the present inven- - 
tion; 

FIG. 2A is a schematic diagram for illustrating dis- 
torted flows of a fluid flowing through a pipe; 

FIG. 2B is a view showing graphically a frequency $0 
distribution of distorted flows; 

FIG. 3 is a flow chart for illustrating a typical exam- 
ple of procedure for calculating a mean flow rate ac- 
cording to the teaching of the present invention; 

FIG. 4 is a view showing an interna] structure of a 65 
memory used for fuzzy inference of the mean flow rate; 

FIG. 5 is a view showing a structure of a flow veloc- 
ity level membership function table; 



FIGS. 6A to 6C are views showing frequency distri- 
butions at different flow velocity levels for illustrating 
the contents of the table shown in FIG. 5; 

FIG. 7 is a view showing a structure of a distorted 
flow model membership function table; 

FIGS. 8A to 8C show three exemplary configura- 
tions of pipe sections through which a fluid under mea- 
surement flows; 

FIGS. 9A to 9J are views showing frequency distri- 
butions of flow velocities in distorted flow model pat- 
terns in dependence on particular shapes of pipes; 

FIG. 10 shows a structure of a mean flow rate calcu- 
lating expression table T3; 

FIG. 11 is a flow chart for illustrating an inference 
procedure executed by a mean flow rate inference part 
of the system; 

FIG. 12 shows a structure of the flow velocity level 
membership function table for illustrating operation of 
the flow rate measuring system according to a second 
embodiment of the invention; 

FIG. 13 is a view showing a structure of the memory 
employed in the flow rate measuring system according 
to a third embodiment of the invention; 

FIG. 14 is a flow chart for illustrating the inference 
procedure executed by a mean flow rate fuzzy inference 
part of the system according to the second embodiment 
of the invention; 

FIG. 15 is a view showing another example of instal- 
lation of sensors within a pipe; 

FIG. 16 is a view showing a structure of a fuzzy 
arithmetic unit according to a fourth embodiment of the 
invention; and 

FIG. 17 is a view showing a general arrangement of 
an air suction system in an internal combustion engine 
to which an embodiment of the flow rate measuring 
system according to the invention is applied. 

DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 

Before entering into detailed description of preferred 
embodiments of the flow rate measuring method and 
system according to the present invention, elucidation 
will first be made of the principle or concept underlying 
a typical example of the mean flow rate calculation 
method taught by the invention. 

FIG. 2A is a schematic diagram for illustrating gener- 
ation of distorted flows in a fluid flowing through a 
pipe. In FIG. 2A, a reference character A denotes gen- 
erally a pipe, tube or conduit (hereinafter generally 
referred to as the pipe) through and within which a fluid 
under measurement flows, and B denotes distorted 
flows of the fluid brought about within the pipe A under 
the influence of geometrical factors such as shape of a 
pipe section located upstream. Further, a reference 
character C denotes generally a flowing direction of the 
fluid. Now, let's consider a flow velocity measurement 
of a fluid flowing through the pipe line A. So long as the 
mean flow rate of the fluid within the pipe line A is 
constant, the ideal conditions for the measurement of 
flow velocities at measuring points Pi, P2 and Pa within 
the pipe A should be such that flow velocity values ui, 
U2 and U3 measured at the respective measuring points 
Pi, P2 and P3 remain always constant with variance 
among the frequency distributions of the flow velocity 
measurement values ui, U2and U3 being zero. However, 
in practical flow velocity measurements, the fluid flow 
is inevitably accompanied with deviations and infinites- 
imal deviations under the influence of disturbances such 
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as vortexes, fluidal separation and/or the like phenom- Referring to FIG. 3, execution of the procedure starts 

ena taking place upstream of the measuring points as from a step SO. In a step SI, detected or measured flow 

well as under the influence of the geometrical configu- velocity values uy U2, . . . , u^/are fetched from flow 

ration of a pipe section located upstream. As a conse- sensors disposed at the measuring points Pim . . . , 

quence, the actual flow velocity (vi, v 2 , v 3 ) exhibits 5 p*,, respectively. Subsequently, in a step S2, a flow 

such frequency distributions which differ from_ one velocity mean (average) value u of ui, U2, . . . , uat is 

another in respect to the temporal mean values vj, v 2 determined on the basis of the flow velocity measure- 

and v 3 and variance, as illustrated in FIG. 2B. In this m ent values ui.nj um in accordance with 

conjunction, it is noted that differences among the tem- 
poral mean values of the flow velocities vj, V2 and V3 10 M f n 
characteristically represent distribution of distorted u = 2 u/Af 
flows, while variance at the individual measuring points . , . 
characterizes the deviation or fluctuation of the flow Ncxt » in 8 ste P S3 > a n ° w veJocit y measurement value 
velocity brought about by the disturbances or turbu- range corresponding to a flow velocity classification 
lences such as vortexes, fluidal separation and/or the 15 ran * e to which the mean velocity of the fluid flow 
like phenomena. When the mean flow rate within the thc P'P* A Wongs and a distorted flow model 
pipe is to be determined on the basis of a plurality of the membership function table bearing a strongest or best 
flow velocity measurement values m, u 2 and u 3 exhibit- correspondence relation to the mean flow velocity are 
ing the distributions described above, it may be con- determined on the basis of the abovementioned mean 
ceived to adopt a method of determining a simple mean 20 flow velocity value u in accordance with thc following 
value of the measured flow velocities ui, U2 and U3 or a expression (2). 

method of determining a weighted flow rate mean value mutVahTHENl^lu (2) 
(wivi+w 2 v 2 +W3V3) by weighting the measured flow 

velocities ui, U2 and U3 with coefficients wi, w 2 and W3, where 

respectively, which can be derived from Reynolds num- 25 * represents a constant in a range of 1, ...» K, 

ber (Re) representing an approximate value of the mean U " ^presents a flow velocity classification range to 

flow rate by taking into account the geometrical config- , which the mean fl ? w vdocit y value u belongs, and 

uration of the flow path and other factors. However, 1 represents a specified or particular range of flow 

the simple mean method is disadvantageous in that flue- velocity measurement values, 

tuations or small excursions will inevitably be included 30 ln order t0 rflow thc distributions of distorted flows 

admixelly in the mean value as determined, because no 10 be handled fuzz V variables, the flow velocity mea- 

consideration is paid to variance of the flow velocity surement values ui, u 2 uj*are divided by a maxi- 

due to the infinitesimal variation of fluctuation of the mum flow velocity value u max given by 
flow velocity. On the other hand, the weighted mean 

method will encounter difficulty in determining defi- 35 <w=m«(«,. uj, . . . . <w) (3) 

nitely the distribution of distorted flows because of the 4 , . u , , - . . 

variance of the flow velocity, which in turn means that t0 * crc * y °* ai! ! normalaed flow velocity measure- 

it becomes difficult to determine evaluation coefficients <? »• U2 ' " ' " u "] ^P 5 ?* and S5 >* A 

related to the distribuUon profile or pattern of the dis- v Subsequently, in a step S6, it is determined on the 
torted flows which are to be taken into consideration in 40 b f «? of * ^ inference rules mentioned below which 

determinating the applicable or pertinent weights w lt ° f the distorted flow model patterns the normalized 

w 2 and W3 flow velocity measurement values ui, 62, . . . , um bears 

At this juncture, there may be mentioned two main the bat rimihrity of fitness (membership grade). TheN 

factors which make it difficult to calculate the mean ^ inference rules are prepared m correspondence to 
flow velocity They are* 45 SP^" 10 » ow velocity levels. More specifically, N dis- 

(1) fuzziness of the mean flow rate brought about by torte * flow t model P atterns " «> m bined with plural 
variations in the mean flow rate level within the pipe, < M > flow velocitv measurement values in the form of a 
and table (refer to, for example, FIG. 7), whereon member- 

(2) fuzziness of the mean flow rate due to variance in or fitness ^ ades (membership function val- 
the flow velocity measurement values due to the dis- 50 u«) are allocated to the abovementioned combinations, 
torted flow distribution and the disturbances such as respectively. 

mentioned above. Expressing mathematically, 

Accordingly, by resorting to the use of a fuzzy infer- 
ence method for determining the mean flow velocity ^ wE q, v 
level and the profile or pattern of the distorted flow 55. 
distribution to thereby take into account the fuzziness in n' ~BQ- Vah 
determination of the weighting coefficients mentioned IF 
above, it is believed to make steady (or stabilize) the 

calculation of the average or mean flow rate, as will be *m -EQ- Vtuk) 
described below in detail 60 

Representing in general by M the number of points at THEN 

which the fluid flow velocity is to be measured, descrip- *~M*nAvi, (4) 
tion will be made on a method of determining the mean 

flow velocity level by resorting to a fuzzy inference where 1= 1, 2, . . . , N, and Jt= 1, 2, . . . , K. 
method. FIG. 3 is a flow chart illustrating a typical 65 In the above expression (4), vi, V2, . . . , vj^are fuzzy 

procedure for calculating a mean fluid flow rate by variables representing typically the flow velocity mea- 

resorting to a fuzzy inference method according to an surement values obtained at the sensor locations or 

embodiment of the invention. measuring points Pi, P2, . . . , Pj* within the pipe, which 
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are normalized by the method described previously, and rate calculating function fa are averaged by weighting 

V/U, Vnkf • * * t V/AT& are fuzzy variables representing the values with the similarity or fitness grade (member- 

the fuzziness of the flow velocity distributions at the ship grade) 4>/jt for the i-th distorted flow model pattern 

abovementioned sensor locations Pi, Pi, . . . t Pat, re- (D m /) of the flow velocity measurement values (ui, U2, . 

spectively, which are derived from the frequency distri- 5 . . , ua/) at the k-th flow velocity level U*.' 

butions of flow velocities at the sensor locations at a Next, on the basis of the mean flow rate tin* deter- 

given level of the mean flow velocity within the pipe mined as described above, a mean value uo is determined 

and normalized through division by a maximum distri- m accordance with the undermentioned expression (7) 

bution frequency. Further, f f * represents a functional t>y using as the weight the membership grade {«for the 

expression for calculating the mean flow rate within the to t-th flow level determined previously (a step S8), 

pipe on the basis of the flow velocity measurement ^hereon the calculation processing comes to an end 

values (ui, U2» . . . ( ua/), wherein the functional expres- ( stcp 
sion may include empirical rules or expressions. More 

specifically, fa represents an expression for calculating K ^ 

the mean flow rate corresponding to the maximum is i ^ k mk 

distribution frequency in the distorted flow distributions « 0 = 

as the i-th fuzzy inference rule, 2 



In the expression (4), "IF' part is referred to as the 

condition part with "THEN" part as the conclusion . . 

part. Notation «.EQ." does not represent equality of 2 Q determination of the mean flow rate in accordance 

values but represents an arithmetic operation for check. with the above expression (7), the weights {* for the 

ing the membership grade of the fuzzy variable v to the individual flow velocity levels and the weights for the 

fuzzy variable V, which will be described below in membership grades to the individual distorted flow 

detail, model patterns are synthetically taken into consider- 

It is now assumed that the flow velocity measurement 2 $ ation * which m turn means that infcrcncc of the mcan 
value ui detected at the measuring point P| is outputted flow rates for intermediate distorted flow distributions 
at a given time point In that case, the membership which are not contained in the distorted flow model 
function indicating to what extent the flow velocity pattern classification table can be realized as well. Fur- 
measurement value ui contributes to validation of the ther, by virtue of the nature of the fuzzy inferences 
condition part, i.e. the grade of equality or fitness of the 30 given by the expressions (4) and (5), unwanted influence 
value ui to the fuzzy variable V,a* is represented by \iu t0 tne inferred value of the mean flow rate can be initi- 
al), where ui represents the normalized flow velocity gated even when any one of the flow velocity measure- 
measurement value at the measuring point Pi. Further, ment values exhibits significant fluctuations due to the 
the membership function of the condition part for the disturbances, because the membership grade (weight) 
flow velocity measurement value (ui, U2, . . . , ujy) is 35 is decreased correspondingly, 
defined as follows: Now, the present invention will be described in detail 

and concrete in conjunction with preferred or exem- 

(5) plary embodiments thereof. 

a. = M :jj u \ FIG. 1 is a block diagram showing schematically a 
* Ul 40 general arrangement of a flow rate measuring system 

according to a first embodiment of the invention. In this 
where a symbol * W represents an operand indicating figure, a reference numeral 1 denotes a fuzzy arithmetic 
that the membership function assumes a minimum (operation) unit which incorporates therein a processor 
value. On these conditions, it is possible to determine unit 11 and a memory 12, wherein the processor unit 11 
with the expression (4) to what extent the flow velocity 45 and the memory 12 are interconnected through an inter- 
distribution of the measurement values (ui, U2, . . . , mm) face 30 so as to allow data transfer therebetween. The 
bears "similarity or fitness" to the distorted flow model memory 12 may be constituted by a random access 
pattern classified as the i-th distorted flow model pat- memory or a read-only memory. On the other hand, the 
tern, when judged on the basis of M flow velocity mea- processor unit 11 includes a central processing unit 

surement values ui mm* wherein given by the 5Q (CPU) 110, a memory 112, an input/output (I/O) circuit 

expression (5) indicates the extent or grade of "similar- 114 and a bus 116 and is in charge of overall control of 

ity or fitness", i.e. the membership grade. In this con- the whole system. Further, the processor unit 11 is 

nection, a minimum value of the term /i^a appearing in provided with a data bus 1110 for outputting externally 

the expression (5) means that the similarity or member- the results of calculations. Additionally, there are pro- 
ship grade is evaluated on the severer conditions. 55 vided for the processor unit 11 input terminals 1111, 

Through similar procedure, the membership grades 4>/* ni2, . . . , 111M which are connected to plural (M in 

are determined for N inference rules, respectively, the illustrated embodiment) sensors 51, 52 5M, 

whereon the mean flow rate within the pipe is calcu- respectively, by way of an analogue-to-digital (A/D) 

lated in accordance with the following expression: converter 2, a noise eliminating filter 3 and respective 

60 detection circuits 41, 42, .... 4M so that the processor 

N (6) unit 11 can fetch data from the sensors 51, 52, ... , 5M, 

/ x i Wnfri. iw) respectively. The sensors 51, 52, . . . , 5M are disposed 

~ n within a pipe 100 through which a fluid flows at the 

i» 1 locations appropriate for detecting distorted flows gen- 

65 erally denoted by a reference numeral 101. In this con- 
where k = 1, 2, . . . , K. junction, it should be mentioned that the sensors 51, 52, 
The above expression (6) equivalently means that . . . , 5M may be disposed on a same plane as shown in 
values determined in accordance with the mean flow FIG. 1 or on mutually different planes as illustrated in 
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FIG. 15 on the condition that these sensors are posi- 
tioned at different heights h. 

The sensors 51, 52, . . . , 5M are constituted by an 
element capable of detecting a change in a physical 
quantity (hereinafter referred to as the physical change) 5 
brought about by motion or kinematic behavior of the 
fluid in which the sensors axe immersed and may each 
be constituted by a pressure sensor, a thermistor; a hot 
wire anemometer or the like. The detection circuits 41, 
42, . . . , 4M serve for function to convert the physical 10 
changes detected by the sensors 51, 52, .... 5M into 
electric signals, respectively. To this end, detection 
circuit may be constituted by a commercially available 
unit selected to conform with the characteristics of the 
sensor. It should further be mentioned that the sensors IS 
51, 52, . . . , 5M may be of such type as to operate on 
mutually different principles with the detection circuits 
41, 42, ... , 4M being then adapted to operate on corre- 
spondingly different principles. 



10 



FIGS. 6A to 6C are views for illustrating, by way of 
example, the contents of the flow velocity membership 
function table shown in FIG. 5. In these figures, there 
are graphically illustrated the flow velocity member- 
ship functions in three different flow velocity ranges 
Urfi, U<ft and Uj3 into which the mean values u of the 
flow velocity measurement values uj, . . . , u m derived 
from the outputs of the sensors 51, 52, . . . , 5M are 
classified (L thus being equal tojhree). They are a low 
flow velocity range U</i (u<fl=u^uji), a medium flow 
velocity range &d\ <u^u<a) and a high flow veloc- 
ity range Urf3(ud2<uSurfj). Referring to FIGS. 6A, 6B 
and 6C taken along the abscissae are three flow veloc- 
ity levels Ui, U2 and U3, respectively, into which the 
mean flow velocity u is classified (with k=3), i.e. a low 
flow velocity level UiJUag^uSuai), a medium flow 
velocity level U2 (u^i < u ^ u^) and a high flow velocity 
level U3 (u a 2<u=u«3), while taken along the ordinates 
the membership values (fitness or membership grades) 



The signals detected by the detection circuits 41, 42, 20 {jt/of the flow velocity levels Ui, U2 and U3 in the flow 



4M are inputted to the noise, eliminating filter 3 
through which high frequency noise components such 
as electric noise, thermal noise and others are elimi- 
nated. Output signals from the filter 3 are converted 
into digital signals by the A/D converter 2 to be subse- 
quently supplied to the fuzzy arithmetic unit 1. 

Through the signal transmission path mentioned 
above, the flow signals detected by the sensors 51, 52, . 
. . , 5M are applied to the input terminals 1111, 1112, 



velocity ranges U<ji, U<a and Ud3» respectively. It will 
be appreciated that the flow velocity level membership 
function table shown in FIG. 5 contains the member- 
ship grades (fitness grades) {wof the flow velocity lev- 
25 els shown in FIGS. 6A, 6B and 6C in the flow velocity 
ranges U<n to Vdt (L=3 in this case), respectively. In 
this manner, when the flow velocity mean value u be- 
longs to or falls within the low flow velocity range U</i, 
the membership grades {n» C21 and {31 of the flow ve- 



, 111M as the signals which can discriminatively be 30 locity levels Ui, U2 and U3 can be determined to be, for 



identified from one another. 

Implemented in the memory 12 are a mean (average) 
flow rate inference part incorporating inference proce- 
dure and tables which are required for the processor 11 
to perform fuzzy inference of the mean flow rate within 
the pipe A on the basis of the signals applied to the input 
terminals 1111, 1112, . . . , 111M. FIG. 4 is a view show- 
ing an internal structure of the memory 12. Referring to 
the FIG. 4, there are prepared in the memory 12 dis- 



example, 1.0, 0.5 and 0.2, respectively, from the table 
content shown in FIG. 6A. Parenthetically, data illus- 
trated in FIGS. 6A to 6C are only for the illustrative 
purpose. It will be understood that L and K may assume 
35 given positive integers, respectively. 

FIG. 7 shows a typical structure of the distorted flow 
model membership function tables T2. Referring to this 
figure, the distorted flow model membership function 
tables Tl includes K tables T21 to T2K which corre- 



torted flow model membership function tables T2 (see 40 spond to the flow velocity ranges Ui to U*, respec- 



F1G. 7) which are required for allowing the profiles of 
distorted flow distributions to be handled as the fuzzy 
variables, a flow velocity level membership function 
table Tl (FIG. 5) for inferring approximately a flow 
velocity range of interest, mean flow rate calculating 45 
expression tables T3 (FIG. 10) for calculating the mean 
flow rate in correspondence to the distorted flow model 
pattern determined through the inference (FIG. 10) and 
a mean flow rate inference memory part Al for storing 
fuzzy inference procedures. In more concrete, the mean 50 
flow rate inference memory part Al stores therein the 
conditional decision statements and the inference ex- 
pressions (1) to (7) described hereinbefore. 

FIG. 5 shows a structure of the flow velocity mem- 
bership function table Tl; Referring to the figure, sym- 55 
bols Urfi, Uj2, . . . , VdL entered in the topmost row of 
the table represent, respectively, specified or particular 
flow velocity ranges resulting from classification of the 
mean values of the measured flow velocities. On the 



tively, which are contained in the flow velocity level 
membership function table Tl, as can be seen in FIG. 4. 
Shown in FIG. 7 is one of these distorted flow model 
membership function tables T21 to T2K, e.g. the table 
T2k for the flow velocity level U* (k is an integer in a 
range given by 1 <k<K). In the table shown in FIG. 7, 
symbols ui, . . . 1 in the topmost row represent 
normalized values of the flow velocity measurement 
values detected by the sensors 51, 52, ... , 5M, respec- 
tively. On the other hand, symbols D m i, D m 2. . . . , 
DmN in the leftmost column represent distorted flow 
model patterns used for classification of the distorted 
flow distribution profiles. Further, symbols (where 
is 1, ... v N and j = l, . . . , M) represent the membership 
function values used for determining the membership 
grades of M flow velocity measurement values to the 
abovementioned distorted flow model patterns, respec- 
tively. As a concrete example of the method for prepar- 
ing the membership grades u,/y*, it is conceivable to 



other hand, symbols Ui, U2, . . . , in the leftmost 60 describe them in the form of interpolation functions by 



column of the table represent flow velocity levels for 
establishing correspondence of the abovementioned 
flow velocity ranges to the distorted flow model mem- 
bership function tables T2. Further, {jt/(k=l, . . . , K 
and 1= 1, . . . , L) represent membership function values 
for establishing relations between the flow velocity 
levels and the distorted flow model membership func- 
tion tables T2. 



making use of the frequency distribution of the flow 
velocity measurement values experimentally deter- 
mined previously or utilize approximate function for 
approximating the mean values and variances of the 
65 flow velocity distributions as the representative param- 
eters by using a Gaussian curve and storing the func- 
tions as subroutines. In this case, it is necessary to nor- 
malize the membership function values by dividing the 
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flow velocity frequency distributions with a maximum 
frequency thereof in conformance with the definition of 
the membership function. 

Now, description will be made of the distorted flow 
model membership function for preparing the distorted 5 
flow model membership function table T2. FIGS. 8A to 
8C show three typical configurations of the flow path 
which are located upstream of the sensors and which 
provide causes for the generation of plural (three) basic 
distorted flow models. In this case, it is assumed, by way 10 
of example, that three flow sensors 5A, 5B and 5C are 
installed. 

FIGS. 9A to 91 are views for illustrating membership 
functions of the normalized values u^, u^and ucof flow 
velocity measurement values va, us and ucoutputted by 1* 
the individual sensors 5A, 5B and 5C in response to a 
distorted flow model D m \ generated by a bent pipe 
having a configuration shown in FIG. 8A. Parentheti- 
cally, the distorted flow model is determined in depen- 
dence on the type of fluid of interest, the geometrical 20 
shape of the pipe or conduit, diameter thereof and so 
forth. Of the figures, FIGS. 9A to 9C illustrate member- 
ship functions at a low flow velocity level Ui resulting 
from classification of the flow velocity into, for exam- 
ple, three levels, FIGS. 9D to 9F illustrate membership 25 
functions at an intermediate flow velocity level U2 and 
FIGS. 9G to 91 illustrate membership functions at a 
high flow velocity level U3. As will be seen in these 
figures, when the mean flow rate u corresponds to the 
low flow velocity level Uj, the membership value (i.e. 30 
fitness or membership grade) \l\a\ corresponding to the 
normalized value uaooT the flow velocity measurement 
value ma detected by the sensor 5A is determined to be 
fta from the graph shown in FIG. 9A. In this way, the 
table T2K includes sub- tables containing the member- 35 
ship function values of the normalized flow velocity 
measurement values ui, . . . , um available from the 
outputs of the sensors 51, 52, , . . , 5M for the distorted 

flow models Dm, D m 2 D m ^for a flow velocity 

model U*. The other tables T21 etc. can be prepared in 40 
the similar manner. 

It should here be mentioned that in place of the dis- 
crete membership function values shown in FIGS. 9A 
to 91, linear interpolation functions such as one shown 
in FIG, 9J may be employed. 45 

FIG. 10 shows a structure of the mean flow rate 
calculating expression table T3. This table includes K 
tables T31 to T3k which correspond to the flow veloc- 
ity levels Ui to U* in the table Tl, as can be seen in FIG. 
4. Of these tables, FIG. 10 shows, by way of example, a 50 
table T3k (where k represents an integer in a range 
given by 1 <k<K) for the flow velocity level U*. Each 
of the flow rate calculating expression tables T31 to T3k 
contains the mean flow rate calculating expressions f,-* 
(i= 1, , . . , N and k= 1, . • . , K) for determining the mean 55 
flow rate from the flow velocity measurement values 
(ui, U2, ; . . , u fl ) for the distorted flow model patterns 

Dmi D m jvat the flow velocity levels Ui, . . . , Ujr, 

respectively. The expression f,-* may be stored in the 
form of interpolation function by using a polynomial of 60 
the flow velocity measurement values (ui, U2, . . . , ujis) 
for the basic distorted flow models at the flow velocity 
levels, respectively, so that the flow rate having the 
maximum distribution frequency can be calculated. 

The expression f,jt may be given, for example, as foi- & 
lows: 



where a/% and b/jt represent coefficients determined for 
every distorted flow model in the tables T31 to T3K. 

Next, referring to FIG. 4, description will be made on 
mutual relations between the flow velocity membership 
function table Tl, the distorted model membership 
function tables T2 and the mean flow rate calculating 
expression tables T3 as well as a data structurizing 
method. The distorted flow model membership func- 
tion table T2 is divided or classified into plural (K) 
function tables T21, T22, . . . , T2K. Similarly, the mean 
flow rate calculating expression table T3 is classified 
into plural (K) mean flow rate calculating expression 
tables T31 ( T32, . . . f T3K. These K function tables T2 
and K expression tables T3 correspond to the flow 
velocity levels Ui, U2, . . . , resulting from division 
of the flow level by K, respectively. In FIG. 4, lines CI, 
C2, . . , , CK represent the abovementioned correspon- 
dence relations, respectively, In the flow velocity level 
membership function table Tl, correspondences are 
established between the measured flow velocity ranges 

U j, U2 Ui and the function tables T21, T22 

T2K and the calculating expression tables T31, T32, . . 
. , T3K through the medium of the weights {*/. Further, 
each of the distorted model membership function tables 
T21, T22, . . . , T2K is divided into plural (N) distorted 
flow model patterns, wherein each of the distorted flow 
model patterns include M membership functions so as to 
be referred to for each of the measurement flow veloc- 
ity values derived from the plural (M) sensor locations. 
Let's consider, for example, the table T2K shown in 
FIG. 4. This table T2K includes plural or N tables 
T2K1, T2K2, .... T2KN, wherein the table T2KN 
contains plural or M distorted flow model membership 
functions pNmk (where m= 1, . . . , M) with correspon- 
dences being established between each of the mean flow 
rate calculating expression tables T31, T32, . . . , T3K 
and each of the plural (N) distorted flow model pat- 
terns. The table T3K, for example, includes plural or N 
mean flow rate calculating expressions f/* (where i= 1, . 

. N), 

The mean flow rate inference part Al is designed to 
infer the mean flow rate of a fluid flowing through a 
pipe in accordance with inference procedure incorpo- 
rated therein by referring to the flow velocity level 
membership function table Tl, the distorted flow model 
membership function tables T2 and the mean flow rate 
calculating expression tables T3 structured as described 
above. In FIG. 4, an arrow RT1 indicates that the mean 
flow rate inference part Al refers to the flow velocity 
level membership function table Tl, and arrows Rl, R2, 
. . . , RK indicate that the mean flow rate inference part 
Al refers to the distorted flow model membership func- 
tion tables T2 and the mean flow rate calculating ex- 
pression tables T3. 

FIG. 11 is a functional block diagram of a calculating 
or computing system for illustrating a procedure 
through which the mean flow rate inference part Al 
infers the mean flow velocity as well as interdepen- 
dence relations among variables. In the figure, arrows 
indicate the interdependence relations among the vari- 
ables. A notation M A1-+TT* represents that the mean 
flow rate inference part Al refers to the flow velocity 
level membership function table Tl. Similarly, symbols 
"AWT2" and M A1-»T3" represent reference of the 
tables T2 and T3 by the inference pan Al, respectively. 
Since description has already been made on the method 
of determining the mean flow rate, the following de- 
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scription wfll be directed to the arithmetic operation ( U*by repeating the steps F7 to F9 and calculates the 

procedure of the mean flow rate inference part Al with mean value uo from the mean values uo* weighted with 

emphasis being put on the manner in which the mean the membership grades (tin accordance with the infer- 

flow rate inference part Al refers to the flow velocity ence equation (7). the mean value uo thus determined is 

level membership function table Tl, the distorted flow 5 supplied to the processor unit 11. 

model membership function tobies T2 and the mean At a step Fll, the arithmetic operation of the infer- 

flow rate calculating expression tables T3. For simplifi- ence part All comes to an end; 

cation of the description, the mean flow rate inference Parenthetically, it should be added that in the first 

part Al is simply termed the inference part Al, while embodiment of the invention described above, the num- 

the flow velocity level membership function table Tl, 10 ber of the distorted flow model patterns may be one (i.e. 

the distorted flow model membership function table T2 jsj a- \y 

and the mean flow rate calculating expression table T3 Next| a embodiment of the invention will be 

are simply called as the tables Tl, T2 and T3, respec- described 

tively. In FIG. 11, symbols F0 to Fll denote step identi- In conjunction with software internal structure 

■ 1C F!' - . m , . „ 15 shown in FIG, 4, it should be mentioned that the mean 

Now referring to FIG. 1, the inference part Al starts flow rate can be inferred by using only the distorted 

operation at a step F0. flow model membership function tables T2 while omit- 

At a step Fl, the inference part Al fetches the flow ti refercnce to the flow vclodty lcvc , membership 

velocity me^urement values^, . , um) fiom the function uWc T1 without nccd for modifyijlg fa M £ 

sensors 51, . . , SM^through the I/O circuit 114 and the 20 ^ expressions of the inference part Al shown in FIG. 

mput termmals 1111 111M (corresponding to the 4> ^ es £ ncc> Ms oomfKmds to that L is set to equal K 

step in mu. 3J. _ A « jA , in the flow velocity level membership function table Tl 

At a step F2 the inference part Al determines the shown in pjQ 5 , herein 
mean value u of the flow velocity measurement values 

(ui, U2, . . . , mm) in accordance with the expression (1) 25 

(refer to the step S2 shown in FIG. 3). &r « l <* - J) (8) 

At a step F3, the inference part Al refers to the table i <* ?t /) 
Tl to thereby determine K membership grades (*/ cor- 
responding to the measured flow rate range (where FIG. 12 shows a corresponding structure of the flow 
1 represents one of 1 to L and k = 1, .... K) to which the 30 velocity level membership function table Tl. In prac- 
mean value u belongs (corresponding to the step S3 tice, however, upon execution of the mean flow rate 
shown in FIG. 3). inference by the processor unit 11, it is required to re- 

At a step F4, the inference part Al determines the place the conditional decision statement given by the 

measured flow velocity maximum value u max (ui, u 2 , . . expression (2) contained in the mean flow rate inference 

. , ujtf) (corresponding to the step S4 in FIG. 3). 35 part A l by the following statement: 

At a step F5, the inference part Al calculates the ip (5eU dk ) THEN | execute k-th one of inference 

normalized flow velocity measurement values ui, . . . , expressions (4)| 

uji/m accordance with u;«u/Ai«„ (where j=l, . . . , M) whereon the fuzzy inference is performed on the dis- 

(corresponding to the step S5 in FIG. 3). torted now distribution pattern. At this juncture, Uj* 

At a step F6, the inference part Al refers to the table 40 represents the k-th measured flow velocity range to 

T2 to determine the membership function values /itf* on w hich the mean flow rate u in one of the measured flow 

the basis of the values ui, G 2) . . . , uji/for every distorted velocity ranges to \J<tk resulting from the division 

flow model patterns D m i, D m2 , . . , D mN at the flow by K belongs, as described previously. Accordingly, in 

velocity levels^, U 2 , . . . , U*(where i= I N, j= 1, the subsequent execution of the inference, there will be 

. . . , M and k— 1, , . . , K). 45 selected only the distorted flow model membership 

At a step F7, the inference part Al calculates the function table T2k and the mean flow rate calculating 

membership grades <f>/*for every distorted flow model expression table T3k. However, since (when 

patterns D mU D m2 , . . . , D mA 'at a flow velocity level k = i) an d { w== o (when k=#), the steps F3 and F10 

Vk in accordance with shown in FIG. 11 are rendered unnecessary with the 

SO result that uojt determined in the step F9 represents uo. 

M According to the method described above, inference 

fak ~ nijtAvj) can be executed at an increased speed because it is suffi- 

7«1 cient to execute only one of the fuzzy inferences for K 

flow velocity levels, although stability or consistency of 

where i — 1, . . . ( N, and j= 1, . . . , M (refer to the step 55 the inferred values for the flow velocity levels may be 

S6 in FIG. 3). degraded more or less. 

At a step F8, the inference part Al refers to the. table With the second embodiment of the invention de- 
T3 to determine the mean flow rate calculating expres- scribed above, speeding-up of the inference is intended 

sions f/jt for every distorted flow model pattern D w j, by using the tables Tl to T3 and the inference expres- 

Dm2» * • • . Dmtf at a given flow velocity level U* (where 60 sions of the first embodiment as they are. It will be 

i=l, . . . , N, and k=l K). however be appreciated that the tables and the infer- 

At a step F9, the inference part Al determines the ence expressions stored in the memory 12 may initially 
mean value uo* from the expressions f,* weighted with be so arranged as to speed up the inference execution, as 
the membership grade <f>/* in accordance with the infer- described below in conjunction with a third embodi- 

ence equation (6) (corresponding to the step S7 in FIG. 65 ment of the invention. 

3). . FIG. 13 is a view showing a structure of the memory 
At a step F10, the inference pan Al determines the 12 implemented to speed up the inference in accordance 

mean values uo* for all the flow velocity levels Ui, . . . with a third embodiment of the invention, wherein 
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contents of the tables Tl, T2 and T3 are same as those 
shown in FIG. 4. Further, FIG. 14 is a flow chart for 
illustrating the inference procedure according to the 
third embodiment. Referring to the figure, the inference 
part Al calculates a mean value u of the flow velocity 5 
measurement values ui, uz, . . . , uj/to thereby determine 
the measured flow velocity range Vdk to which the 
mean flow velocity value u belongs, as described here- 
inbefore in conjunction with the second embodiment 
(step F2). Subsequently, by referring to only the tobies 10 
T2k and T3k specified by the flow velocity range Vdk as 
determined, the steps F4 to F9 are executed. The value 
nok determined at the step F9 is outputted as the mean or 
average value uo- 

As will now be understood from the foregoing, the IS 
inference methods according to the embodiments of the 
invention do not reside in determining merely a simple 
mean value of the individual sensor outputs but takes 
into account the characteristics of the distorted flow 
pattern prevailing at the sensor locations upon deter- 20 
mining the mean value of the sensor outputs. By virtue 
of this feature, a plurality of the sensors need not neces- 
sarily be located on a same plane transversal to the axis 
of the pipe or within a predetermined range in the direc- 
tion axial! y or diameterwise of the pipe. In other words, 25 
the individual sensors 51, 52, ... , 5M may be disposed 
at different positions along the flow of a fluid within the 
pipe with heights h which differ from one to another 
sensor, as is illustrated in FIG. 15. Thus, there can be 
increased the degree of freedom of selecting the posi- 30 
toons internally of the pipe where the sensors are to be 
installed. With the disposition shown in FIG. 15, the 
individual sensors are disposed at different positions 
along the fluid flowing direction, whereby a temporal 
mean value of the flow rate can be determined. 35 

In order to speed up the fuzzy inference, it is effective 
to execute the weighted mean determining arithmetic 
operation (fuzzy integration) at the steps F2, F9 and 
F10 in the inference procedure shown in FIG. 11. FIG. 
16 shows a fourth embodiment of the invention imple- 40 
mented to this end. Referring to the figure, a numerical 
integration processor 13 incorporated in the fuzzy arith- 
metic unit 1 for performing a numerical integration is 
provided and connected to the processor unit 11 via the 
data bus 31. 45 

FIG. 17 is a view showing an air suction control 
system of an internal combustion engine to which an 
embodiment of the flow rate measuring system accord- 
ing to the invention is applied. Referring to the figure, a 
suction air flow denoted by a reference symbol SA is 50 
forced to flow through an air filter 503 contained in an 
air cleaner 511, a suction pipe 100, an air flow meter 
102, a throttle valve 103 and an intake manifold 501 to 
be charged into a combustion chamber or chambers of 
the internal combustion engine. 55 

Disposed in a main flow path of the air flow meter 
102 are flow velocity sensors 51, 52, . . . , 5M which are 
connected to detection circuits 41, 42, ... , 4M, respec- 
tively, whose outputs are in turn connected through a 
noise eliminating filter 3 and an analogue-to-digital 60 
(A/D) converter 2 to a fuzzy arithmetic unit 1 which is 
implemented in accordance with one of the first to third 
embodiments of the invention described above and by 
which the mean flow rate of air flowing through the 
suction pipe 100 is calculated by taking into account the 65 
distorted flow components generally denoted by 102. 
The throttle valve 9 disposed at a position downstream 
of the air flow meter 102 is interlocked to an accelerator 



pedal (not shown) for controlling the suction or intake 
air flow. Mounted on the throttle body 103 is an idle 
speed control valve 8 which serves for controlling the 
air flow at the time when the throttle valve 9 is fully 
closed. In FIG. 17, thick solid arrows indicate direc- 
tions of fuel flows. The fuel pumped up from a fuel tank 
505 by a fuel injection pump 506 is injected into within 
the manifold 501 through fuel injectors 507 to be mixed 
with the air flow having passed through the air flow 
meter 103, the resulting air-fuel mixture being charged 
into the engine 500. 

A control unit 510 serves for arithmetically determin- 
ing the amount of fuel injected through the injector and 
the opening degree of the idle speed control valve 8 on 
the basis of the output of the fuzzy arithmetic unit 1, an 
angular signal representing rotation of the throttle valve 
9, an output signal of an oxygen concentration sensor 
508 disposed within an exhaust manifold 511 and an 
output signal of an engine rotation speed sensor 509. In 
accordance with the results of the fuzzy arithmetic 
operation, the injectors 507 and the idle speed control 
valve 8 are correspondingly controlled. 

The air suction pipe system of the internal combus- 
tion engine is often of a much intricate structure because 
of little space available for installation of the piping in 
the engine room, which makes it difficult to provide an 
auxiliary air flow path. Consequently, the air flow 
within the suction pipe contains significant disturbances 
or turbulence, fluidal noise and distorted flow compo- 
nents. By incorporating the flow measurement appara- 
tus according to the invention in this kind of suction 
piping, the mean flow rate can be measured consistently 
with significantly improved accuracy regardless of the 
presence of distorted flows and noise, whereby high- 
precision air-fuel ratio control can be realized. 

As will be apparent from the foregoing, the present 
invention provides advantageous effects mentioned 
below. 

(1) Because of use of fuzzy integration for calculation 
of the intra-pipe mean flow rate, the flow rate value as 
determined is not affected by the output of any particu- 
lar sensor, whereby the mean flow rate value can be 
obtained with high consistency and reliability, being 
protected against influence of local flow disturbances. 

(2) Since the membership grade or fitness grade of a 
distorted flow pattern is determined for each of plural 
classified distorted flow model patterns and since ex- 
tents of contribution of the distorted flow model pat- 
terns to the determination of the mean flow rate by 
using the membership grades (fitness grades) as the 
weights, it is possible to calculate appropriately the 
mean flow rate even for those distorted flow distribu- 
tions which are intermediate the distorted flow model 
patterns. 

(3) The membership functions for the distorted flow 
model patterns can be changed or modified by exchang- 
ing software. As a result, the flow rate measurement can 
flexibly be performed for a variety of distorted flows 
which differ from one another in respect to the distribu- 
tion and the position where the disturbance takes place 
with the sensors being fixed in respect to the shape, 
location of installation, characteristics and other hard- 
ware aspects. 

(4) The flow rate measuring system according to the 
invention can profitably be employed in internal com- 
bustion engines of motor vehicles. In that case, since the 
air flow rate can be measured with a high accuracy even 
in an air piping suction system of an intricate structure, 
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the fuel injection amount can be regulated finely in 
dependence on the measured air flow rate. Thus, an 
air-fuel ratio control assuring a very high accuracy can 
be realized. In general, the geometrical factors such as 
shape of the pipe in the air suction system vary in depen- 
dence on the types of the motor vehicles. In this con- 
junction, the flow rate measuring system according to 
the invention can be applied to any types of the vehicles 
without need for any modifications including that of 
data contained in the memory. 

The teachings of the invention can be applied to flow 
rate measurement of a variety of fluids regardless of 
whether it is liquid or gas. 

I claim: 

1. A flow rate measuring method, in control appara- 
tus for controlling an internal combustion engine, of 
determining an intake flow rate of a fluid flowing 
through a suction pipe of said internal combustion en- 
gine by using a fuzzy inference system which is imple- 
mented on the basis of fuzzy knowledge including fuzzy 
rules each described in the form of an IF part and a 
THEN part, and membership functions defining mean- 
ings of statements described in said IF and THEN parts, 
said control apparatus including speed detecting means 
for detecting rotational speed of said internal combus- 
tion engine, fuel injection means for injecting a fuel into 
the intake air flow within said suction pipe and control 
means for regulating said fuel injection amount of said 
fuel injection means on the basis of weighted flow rate 
means value of said fluid and output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine said method 
comprising steps of: 



10 



IS 



20 



25 



30 



speed of said internal combustion engine, said method 
comprising steps of: 
measuring flow velocities of said fluid at a plurality of 

locations within said pipe; 
preparing a membership function table containing 
membership functions indicating grades of fitness 
of said flow velocities to each of plural distorted 
flow model distributions; 
preparing mean flow rate calculating expressions 
giving flow rates corresponding to said distorted 
flow model distributions, respectively; 
determining membership grades of the measured flow 
velocities to each of said distorted flow model 
distributions by referring to said membership func- 
tion table through the medium of said measured 
flow velocities; and 
calculating said weighted mean of said flow rates 
given by said mean flow rate calculating expres- 
sions by using as weights said membership grades 
to said distorted flow model distributions. 
3. A flow rate measuring method, in control appara- 
tus for controlling an internal combustion engine, of 
determining an intake flow rate of a fluid flowing 
through a suction pipe of said internal combustion en- 
gine by using a fuzzy inference system which is imple- 
mented on the basis of fuzzy knowledge including fuzzy 
rules each described in the form of an IF part and a 
THEN part, and membership functions defining mean- 
ings of statements described in said IF and THEN parts, 
said control apparatus including speed detecting means 
for detecting rotational speed of said internal combus- 
tion engine, fuel injection means for injecting a fuel into 
the intake air flow within said suction pipe, and control 



measuring flow velocities of said fluid at a plurality of 35 means for regulating said fuel injection amount of said 

locations within said suction pipe; fuel injection means on the basis of weighted mean of 

preparing a table of membership functions indicating flow rates of said fluid and output signal of said speed 

grades of fitness of said flow velocities to a dis- detecting means for thereby controlling the rotational 

torted flow model distribution; speed of said internal combustion engine, said method 

preparing a mean flow rate calculating expression comprising steps of: 



giving a flow rate corresponding to said distorted 
flow model distribution; 
determining a membership grade of the measured 
flow velocities to said distorted flow model distri- 
bution by referring to said membership function 45 
table through the medium of said measured flow 
velocities; and 
calculating said weighted flow rate mean value given 
by said mean flow rate calculating expression by 
using as a weight said membership grade to said 
distorted flow model distribution. 
2. A flow rate measuring method, in control appara- 
tus for controlling an internal combustion engine, of 
determining an intake flow rate of a fluid flowing 
through a suction pipe of said internal combustion en- 
gine by using a fuzzy inference system which is imple- 
mented on the basis of fuzzy knowledge including fuzzy 
rules each described in the form of an IF part and a 
THEN part, and membership functions defining mean- 
ings of statements described in said IF and THEN parts, 
said control apparatus including, speed detecting means 
for detecting rotational speed of said internal combus- 
tion engine, fuel injection means for injecting a fuel into 
the intake air flow within said suction pipe, and control 
means for regulating said fuel injection amount of said 65 
fuel injection means on the basis of a weighted means of 
flow rates of said fluid and an outut signal of said speed 
detecting means for thereby controlling the rotational 



50 



55 



60 



measuring flow velocities of said fluid at a plurality of 

locations within said pipe; 
preparing a first membership function table indicating 
grades of fitness of a mean value of said flow veloc- 
ities within said pipe to each of plural flow velocity 
levels; 

preparing second membership function tables giving 
grades of fitness of said flow velocities to a dis- 
torted flow model distribution at each of said flow 
velocity levels; 

preparing mean flow rate calculating expressions 
giving flow rates corresponding to said distorted 
flow model distribution at said flow velocity levels, 
respectively; 

determining a mean value of the flow velocities mea- 
sured at said plurality of locations within said pipe 
to thereby determine on the basis of said mean 
value first membership grades of said mean value to 
said flow levels, respectively, by reference to said 
first membership function table through the me- 
dium of said mean value; 

determining a second membership grade of said flow 
velocity values to said distorted flow model distri- 
bution at each of said flow velocity levels by refer- 
ence to said second membership function tables 
through the medium of said measured flow veloc- 
ity values; and 

calculating said weighted mean of said flow 
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rates given by said mean flow rate calculating expres- 
sions for said flow velocity levels by using as 
weights said first and second membership grades. 

4. A flow rate measuring method, in control appara- 
tus for controlling an internal combustion engine, of 5 
determining an intake flow rate of a fluid flowing 
through a suction pipe of said internal combustion by 
using a fuzzy inference system which is implemented on 
the basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 10 
and membership functions defining meanings of state- 
ments described in said IF and THEN parts, said con- 
trol apparatus including speed detecting means for de- 
tecting rotational speed of said internal combustion 
engine, fuel injection means for injecting a fuel into the 
intake air flow within said suction pipe, and control 
means for regulating said fuel injection amount of said 
fuel injection means on the basis of a weighted mean of 
flow rates of said fluid and an output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine, said method 
comprising steps of: 

measuring flow velocities of said fluid at a plurality of 
locations within said suction pipe; 

preparing a first membership function table contain- 
ing membership functions indicating grades of fit- 
ness of a mean value of said flow velocities within 
said suction pipe to plural flow velocity levels, 
respectively; 3Q 

preparing second membership function tables con- 
taining membership functions indicating grades of 
fitness of said flow velocities to a plurality of dis- 
torted flow model distributions at said flow veloc- 
ity level, respectively; 35 

preparing mean flow rate calculating expressions 
giving flow rates corresponding to said distorted 
flow model distributions at said flow velocity lev* 
els, respectively; 

determining a mean value of the flow velocities mea- 4Q 
sured at said plurality of locations within said suc- 
tion pipe to thereby determine on the basis of said 
mean value first membership grades of said mean 
value to said flow levels, respectively, by reference 
to said first membership function table through the 45 
medium of said mean value; 

determining second membership grades or said flow 
velocity values to said distorted flow model distri- 
butions at said flow velocity levels, respectively, 
by reference to said second membership function 50 
table through the medium of said measured flow 
velocity values; and 

calculating said weighted mean of said flow rates 
given by said mean flow rate calculating expres- 
sions for said flow velocity levels, respectively, by 55 
using as weights said first and second membership 
grades. 

5. A flow rate measuring method according to claim 

4, 

wherein the step of calculating the flow rate within 60 
said pipe includes: 

a step of determining weighted flow rate mean values 
in accordance with said mean flow rate calculating 
expression by using as weights said second mem- 
bership grades as determined; and 65 

a step of deriving a weighted mean of said weighted 
flow rate mean values calculated in accordance 
with said mean flow rate calculating expression for 



said flow velocity levels by using as weights said 
first membership grades determined. 

6. A flow rate measuring system, in control apparatus 
for controlling an internal combustion engine, for deter- 
mining an intake flow rate of a fluid flowing through a 
suction pipe of said internal combustion engine by using 
a fuzzy inference system which is implemented on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and THEN part, and 
membership functions defining meanings of statements 
described in said IF and THEN parts, said control appa- 
ratus including speed detecting means for detecting 
rotational speed of said internal combustion engine, fuel 
injection means for injecting a fuel into the intake air 
flow within said suction pipe, and control means for 
regulating said fuel injection amount of said fuel injec- 
tion means on the basis of an weighted mean of flow 
rates of said fluid and an output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine, said flow rate 
measuring system comprising: 

a plurality of detecting means for measuring flow 
velocities of said fluid at a plurality of locations 
within said suction pipe; 

storage means containing a table of membership func- 
tions giving grades of fitness of said flow velocities 
to a distorted flow model distribution; 

storage means containing a mean flow rate calculat- 
ing expression giving a flow rate corresponding to 
said distorted flow model distribution; 

means for determining a membership grade of the 
measured flow velocities to said distorted flow 
model distribution by referring to said membership 
function table through the medium of said mea- 
sured flow velocity; and 

means for determining a weighted mean value of said 
flow rates given by said mean flow rate calculating 
expression by suing as a weight said membership 
grade to said distorted flow model distribution. 

7. A flow rate measuring system, in control apparatus 
for controlling an internal combustion engine, for deter- 
mining an intake flow rate of a fluid flowing through a 
suction pipe of said internal combustion engine by using 
a fuzzy inference system which is implemented on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 
and membership functions defining meanings of state- 
ments described in said IF and THEN parts, said con- 
trol apparatus including speed detecting means for de- 
tecting rotational speed of said internal combustion 
engine, fuel injection means for injecting a fuel into the 
intake air flow within said suction pipe, and control 
means for regulating said fuel injection amount of said 
fuel injection means on the basis of an weighted mean of 
flow rates of said fluid and an output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine, said flow rate 
measuring system comprising: 

a plurality of detecting means for measuring flow 
velocities of said fluid at a plurality of locations 
within said suction pipe; 

storage means for storing a first membership function 
table containing membership functions indicating 
grades of fitness of a mean value of flow velocities 
within said suction pipe to a plurality of distorted 
flow velocity levels, respectively; 

storage means for storing second membership func- 
tion tables containing second membership func- 
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tions indicating grades of fitness of said flow veloc- 
ities to a plurality of distorted flow model distribu- 
tions at said flow velocity levels, respectively; 

means for determining membership grades of said 
flow velocity values to said distorted flow model 5 
distributions, respectively, by reference to said 
second membership function tables through the 
medium of said measured flow velocity values; and 

means for calculating a weighted mean of said flow 
rates given by said mean flow rate calculating ex- 10 
pressions for said flow velocity levels, respectively, 
by using as weights said membership grades. 

8. A flow rate measuring system, in control apparatus 
for controlling an internal combustion engine, for deter- 
mining an intake flow rate of a fluid flowing through a IS 
suction pipe of said internal combustion engine by using 

a fuzzy inference system which is implemented on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 
and membership functions defining meanings of state- 20 
ments described in said IF and THEN parts, said con- 
trol apparatus including, speed detecting means for 
detecting rotational speed of said internal combustion 
engine, fuel injection means for injecting a fuel into the 
intake air flow within said suction pipe, and control 25 
means for regulating said fuel injection amount of said 
fuel injection means on the basis of a weighted mean of 
flow rates of said fluid and an output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine, said flow rate 30 
measuring system comprising: 
a plurality of detecting means for mesuring flow ve- 
locities of said fluid at a plurality of locations 
within expression by using as weights said second 
membership grades said suction pipe; 35 
storage means containing a first membership function 
table indicating grades of fitness of a mean value of 
said flow velocities within said pipe to a distorted 
flow velocity level; 
storage means containing second membership func- 40 
tion tables indicating grades of fitness of said flow 
velocities to distorted flow model distributions at 
said flow velocity levels, respectively; 
storage means containing a mean flow rate calculat- 
ing expressions giving a flow rate corresponding to 45 
said distorted flow model distribution at each of 
said flow velocity levels; 
means for determining a mean value of the flow ve- 
locities measured at said plurality of locations 
within said suction pipe to thereby determine on 50 
the basis of said mean value a first membership 
grade of said mean value to said flow levels, respec- 
tively, by reference to said first membership Amo- 
tion table through the medium of said mean value; 
means for determining a second membership grade of 55 
said measured flow velocity values to said distorted 
flow model distributions for said flow velocity 
levels, respectively, by reference to said second 
membership function table through the medium of 
said measured flow velocity values; and 60 
means for calculating said weighted mean of said 
flow rates given by said mean flow rate calculating 
expressions for said flow velocity level, respec- 
tively, by using as weights said first and second 
membership grades, for thereby determining the 65 
flow rate of the fluid flowing within said pipe. 

9. A flow rate measuring system, in control apparatus 
for controlling an internal combustion engine, for deter- 
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mining an intake flow rate of a fluid flowing through a 
suction pipe of said internal combustion engine by using 
a fuzzy inference system which is implemented on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 
and membership functions defining meanings of state- 
ments described in said IF and THEN parts, said con- 
trol apparatus including, speed detecting means for 
detecting rotational speed of said internal combustion 
engine, fuel injection means for injecting a fuel into the 
intake air flow within said suction pipe, and control 
means for regulating said fuel injection amount of said 
fuel injection means on the basis of a weighted mean of 
flow rates of said fluid and an output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine, said flow rate 
measuring system comprising: 
a plurality of detecting means for measuring flow 
velocities of said fluid at a plurality of locations 
within said suction pipe; 
storage means containing a first membership function 
table indicating grades of fitness of a mean value of 
said flow velocities within said suction pipe to a 
plurality of distorted flow velocity level, respec- 
tively; 

storage means containing a second membership func- 
tion table indicating grades of fitness of said flow 
velocities to a plurality of distorted flow model 
distributions for said flow velocity levels, respec- 
tively; 

storage means continuing a mean flow rate calculat- 
ing expressions for deriving flow rates correspond- 
ing to said distorted flow model distributions at 
said flow velocity levels, respectively; 
means for determining a mean value of the flow ve- 
locities measured at said plurality of locations 
within said suction pipe to thereby determine on 
the basis of said mean value a first membership 
grades of said mean value to said flow levels, re- 
spectively, by reference to said first membership 
function table through the medium of said mean 
value; 

means for determining a second membership grades 
of said flow velocity values to said distorted flow 
model distributions for said flow velocity levels, 
respectively, by reference to said second member- 
ship function table through the medium of said 
measured flow velocity values; and 
means for calculating said weighted mean of said 
flow rates given by said mean flow rate calculating 
expressions, respectively, by using as weights said 
first and second membership grades, for thereby 
determining the flow rate of the fluid within said 
suction pipe. 

10. A flow rate measuring system according to claim 

wherein said means for calculating the flow rate 

within said pipe includes: 
means for determining weighted flow rate mean val- 
ues in accordance with said mean flow rate calcu- 
lating expression by using as weights said second 
membership grades as determined; and 
means for deriving a weighted mean of said weighted 
flow rate mean values calculated in accordance 
with said mean flow rate calculating expression for 
said flow velocity levels by using as a weight said 
first membership grade as determined. 
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11. A flow rate measuring system according to claim 
9, wherein said plurality of detecting means are dis- 
posed in a same plane extending orthogonally to the 
longitudinal axis of said pipe. 

12. A flow rate measuring system according to claim 5 
9, wherein said plurality of detecting means are dis- 
posed along a direction in which said fluid flows within 
said pipe. 

13. A flow rate measuring method of determining a 
flow rate of a fluid flowing through an intake pipe of an 10 
internal combustion engine of an internal combustion 
engine by using a flow rate measuring system including 

a fuzzy inference system which is implemented on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 15 
and membership functions defining meanings of state- 
ments described in said IF and THEN parts, 
said method comprising steps* performed by said 

flow rate measuring system, of: 
measuring flow velocities of said fluid at a plurality of 20 

locations within said pipe; 
preparing a table of membership functions indicating 
grades of fitness of said flow velocities to a dis- 
torted flow model distribution; 
preparing a mean flow rate calculating expression 25 
giving a flow rate corresponding to said distorted 
flow model distribution; 
determining a membership grade of the measured 
flow velocities to said distorted flow model distri- 
bution by referring to said membership function 30 
table through the medium of said measured flow 
velocities; 

calculating a weighted flow rate mean value given by 
said mean flow rate calculating expression by using 
as a weight said membership grade to said distorted 35 
flow model distribution; and 

regulating fuel injection amount of fuel injection 
amount of fuel injection means of said internal 
combustion engine on the basis of said weighted 
flow rate mean value and a signal indicative of 40 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 

14. A flow rate measuring method of determining a 
flow rate of a fluid flowing through an intake pipe of an 45 
internal combustion engine of an interna] combustion 
engine by using a flow rate measuring system including 

a fuzzy inference system which is implemented on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 50 
and membership functions defining meanings of state- 
ments described in said IF and THEN parts, 
said method comprising steps, performed by said 

flow rate measuring system, of: 
measuring flow velocities of said fluid at a plurality of 55 

locations within said pipe; 
preparing a membership function table containing 
membership functions indicating grades of fitness 
of said flow velocities to each of plural distorted 
flow model distributions; 60 
preparing mean flow rate calculating expressions 
giving flow rates corresponding to said distorted 
flow model distributions, respectively; 
determining membership grades of the measured flow 
velocities to each of said distorted flow model 65 
distributions by referring to said membership func- 
tion table through the medium of said measured 
flow velocities; and 



calculating a weighted mean of said flow rates given 
by said mean flow rate calculating expressions by 
using as weights said membership grades to said 
distorted flow model distributions; and 

regulating fuel injection amount of fuel injection 
amount of fuel injection means of said internal 
combustion engine on the basis of said weighted 
mean of said flow rates and a signal indicative of 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 

15. A flow rate measuring method of detennining a 
flow rate of a fluid flowing through an intake pipe of an 
internal combustion engine of an internal combustion 
engine by using a fuzzy inference system which is im- 
plemented on the basis of fuzzy knowledge including 
fuzzy rules each described in the form of an IF part and 
a THEN part, and membership functions defining 
meanings of statements described in said IF and THEN 
parts, 

said method comprising steps, performed by said 
flow rate measuring system, of: 

measuring flow, velocities of said fluid at a plurality of 
locations within said pipe; 

preparing a first membership function table indicating 
grades of fitness of a mean value of said flow veloc- 
ities within said pipe to each of plural flow velocity 
levels; 

preparing second membership function tables giving 
grades of fitness of said flow velocities to a dis- 
torted flow model distribution at each of said flow 
velocity levels; 

preparing mean flow rate calculating expressions 
giving flow rates corresponding to said distorted 
flow model distribution at said flow velocity levels, 
respectively; 

determining a mean value of the flow velocities mea- 
sured at said plurality of locations within said pipe 
to thereby determine on the basis of said mean 
value first membership grades of said mean value to 
said flow levels, respectively, by reference to said 
first membership function table through the me- 
dium of said mean value; 

determining a second membership grade of said flow 
velocity values to said distorted flow model distri- 
bution at each of said flow velocity levels by refer- 
ence to said second membership function tables 
through the medium of said measured flow veloc- 
ity values; 

calculating a weighted mean of said flow rates given 
by said mean flow rate calculating expressions for 
said flow velocity levels by using as weights said 
first and second membership grades; and 

regulating fuel injection amount of fuel injection 
amount of fuel injection means of said internal 
combustion engine on the basis of said weighted 
mean of said flow rates and a signal indicative of 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 

16. A flow rate measuring method of detennining a 
flow rate of a fluid flowing through an intake pipe of an 
internal combustion engine by using a flow rate measur- 
ing system including fuzzy inference system which is 
implemented on the basis of fuzzy knowledge including 
fuzzy rules each described in the form of an IF part and 
a THEN part, and membership functions defining 
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meanings of statements described in said IF and THEN 
parts, 

said method comprising steps, performed by said 
flow rate measuring system, of: 

measuring flow velocities of said fluid at a plurality of 5 
locations within said pipe; 

preparing a first membership function table contain- 
ing membership functions indicating grades of fit- 
ness of a mean value of said flow velocities within 
said pipe to plural flow velocity levels, respec- 10 
tively; 

preparing second membership function tables con* 
taining membership functions indicating grades of 
fitness of said flow velocities to a plurality of dis- 
torted flow model distributions at said flow veloc- 15 
hy level, respectively; 
preparing mean flow rate calculating expressions 
giving flow rates corresponding to said distorted 
flow model distributions at said flow velocity lev- 
els, respectively; 20 
determining a mean value of the flow velocities mea- 
sured at said plurality of locations within said pipe 
to thereby determine on the basis of said mean 
value first membership grades of said mean value to 
said flow levels, respectively, by reference to said 25 
first membership function table through the me- 
dium of said mean value; 
determining second membership grades or said flow 
velocity values to said distorted flow model distri- 



a plurality of detecting means for measuring flow 
velocities of said fluid at a plurality of locations 
within said pipe; 
storage means containing a table of membership func- 
tions giving grades of fitness of said flow velocities 
to a distorted flow model distribution; 
storage means containing a mean flow rate calculat- 
ing expression giving a flow rate corresponding to 
said distorted flow model distribution; 
means for determining a membership grade of the 
measured flow velocities to said distorted flow 
mode] distribution by referring to said membership 
function table through the medium of said mea- 
sured flow velocity; 
means for determining a weighted mean value of said 
flow rates given by said mean flow rate calculating 
expression by suing as a weight said membership 
grade to said distorted flow model distribution; and 
control means for regulating fuel injection amount of 
fuel injection means of said internal combustion 
engine on the basis of said weighted mean of said 
flow rates mean value and a signal indicative of 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 
19. A flow rate measuring system for determining a 
flow rate of a fluid flowing through an intake pipe of an 
internal combustion engine by using a flow rate measur- 
ing system including fuzzy inference system which is 



butions at said flow velocity levels, respectively, 30 implemented on the basis of fuzzy knowledge including 



by reference to said second membership function 
table through the medium of said measured flow 
• velocity values; 

calculating a weighted mean of said flow rates given 
by said mean flow rate calculating expressions for 35 
said flow velocity levels, respectively, by using as 
weights said first and second membership grades; 
and 

regulating fuel injection amount of fuel injection 
amount of fuel injection means of said internal 40 
combustion engine on the basis of said weighted 
mean of said flow rates and a signal indicative of 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 45 

17. A flow rate measuring method according to claim 
16, 

wherein the step of calculating the flow rate within 
said pipe includes: 

a step of determining weighted flow rate mean value 20 
in accordance with said mean flow rate calculating 
expression by using as weights said second mem- 
bership grades as determined; and 

a step of deriving a weighted mean of said weighted 
flow rate mean values calculated in accordance 55 
with said mean flow rate calculating expression for 
said flow velocity levels by using as weights said 
first membership grades determined. 

18. A flow rate measuring system for determining a 
flow rate of a fluid flowing through an intake pipe of an 60 
internal combustion engine by using a flow rate measur- 
ing system including fuzzy inference system which is 
implemented on the basis of fuzzy knowledge including 
fuzzy rules each described in the form of an IF part and 

a THEN part, and membership functions defining 65 
meanings of statements described in said IF and THEN 
parts, 

said flow rate measuring system comprising: 



fuzzy rules each described in the form of an IF part and 
a THEN part, and membership functions defining 
meanings of statements described in said IF and THEN 
parts, 

said flow rate measuring system comprising: 

a plurality of detecting means for measuring flow 
velocities of said fluid at a plurality of locations 
within said pipe; 

storage means for storing a first membership function 
table containing membership function indicating 
grades of fitness of a mean value of said flow veloc- 
ities within said pipe to a plurality of distorted flow 
velocity levels, respectively; 

storage means for storing second membership func- 
tion tables containing second membership func- 
tions indicating grades of fitness of said flow veloc- 
ities to a plurality of distorted flow model distribu- 
tions at said flow velocity levels, respectively; 

means for deterauning membership grades of said 
flow velocity values to said distorted flow model 
distributions, respectively, by reference to said 
second membership function tables through the 
medium of said measured flow velocity values; 

means for calculating a weighted mean of said flow 
rates given by said mean flow rate calculating ex- 
pressions for said flow velocity levels, respectively, 
by using as weights said membership grades; and 

control means for regulating fuel injection amount of 
fuel injection means of said internal combustion 
engine on the basis of said weighted means of said 
flow rates mean value and a signal indicative of 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 

20. A flow rate measuring system for determining a 
flow rate of a fluid flowing through an intake pipe of an 
internal combustion engine by using a flow rate measur- 
ing system including fuzzy inference system which is 
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implemented on the basis of fuzzy knowledge including 
fuzzy rules each described in the form of an IF part and 
a THEN part, and membership functions defining 
meanings of statements described in said IF and THEN 
parts, 5 
said flow rate measuring system comprising: 
a plurality of detecting means for measuring flow 
velocities of said fluid at a plurality of locations 
within said pipe; 
storage means containing a first membership function 1° 
table indicating grades of fitness of a mean value of 
said flow velocities within said pipe to a distorted 
flow velocity level; 
storage means containing second membership func- 
tion tables indicating grades of fitness of said flow 15 
velocities to distorted flow model distributions at 
said flow velocity levels, respectively; 
storage means containing a mean flow rate calculat- 
ing expressions giving a flow rate corresponding to 
said distorted flow model distribution at each of 
said flow velocity levels; 
means for determining a mean value of the flow ve- 
locities measured at said plurality of locations 
within said pipe to thereby determine on the basis 2J 
of said mean value a first membership grade of said 
mean value to said flow levels, respectively, by 
reference to said first membership function table 
through the medium of said mean value; 
means for determining a second membership grade of 30 
said measured flow velocity values to said distorted 
flow model distribution for said flow velocity lev- 
els, respectively, by reference to said second mem- 
bership function table through the medium of said 
measured flow velocity values; 35 
means for calculating a weighted mean of said flow 
rates given by said mean flow rate calculating ex- 
pressions for said flow velocity level, respectively, 
by using as weights said first and second member- 
ship grades, for thereby determining the flow rate 40 
of the fluid flowing within said pipe; and 
control means for regulating fuel injection amount of 
fuel injection means of said internal combustion 
engine on the basis of said weighted mean of said 
flow rates mean value and a signal indicative of 45 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 
21. A flow rate measuring system for determining a 
flow rate of a fluid flowing through an intake pipe of an 50 
internal combustion engine by using a flow rate measur- 
ing system including fuzzy inference system which is 
implemented on the basis of fuzzy knowledge including 
fuzzy rules each described in the form of an IF part and 
a THEN part, and membership functions defining 55 
meanings of statements described in said IF and THEN 
parts, 

said flow rate measuring system comprising: 

a plurality of detecting means for measuring flow 

velocities of said fluid at a plurality of locations 60 

within said pipe; 



storage means containing a first membership function 
table indicating grades of fitness of mean value of 
said flow velocities within said pipe to a plurality of 
distorted flow velocity level, respectively; 

storage means containing a second membership func- 
tion table indicating grades of fitness of said flow 
velocities to a plurality of distorted flow model 
distributions for said flow velocity levels, respec- 
tively; 

storage means continuing mean flow rate calculating 
expressions for deriving flow rates corresponding 
to said distorted flow model distributions at said 
flow velocity levels, respectively; 

means for determining a mean value of the flow ve- 
locities measured at said plurality of locations 
within said pipe to thereby determine on the basis 
of said mean value a first membership grades of 
said mean value to said flow levels, respectively, by 
reference to said first membership function table 
through the medium of said mean value; 

means for determining a second membership grades 
of said flow velocity values to said distorted flow 
model distributions for said flow velocity levels, 
respectively, by reference to said second member- 
ship function table through the medium of said 
measured flow velocity values; 

means for calculating a weighted mean of said flow 
rates given by said mean flow rate calculating ex- 
pressions, respectively, by using as weights said 
first and second membership grades, for thereby 
determining the flow rate of the fluid within said 
pipe; and 

control means for regulating fuel injection amount of 
fuel injection means of said internal combustion 
engine on the basis of said weighted mean of said 
flow rates mean value and a signal indicative of 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 

22. A flow rate measuring system according to claim 
21, 

wherein said means for calculating the flow rate 
within said pipe includes: 

means for determining weighted flow rate mean val- 
ues in accordance with said mean flow rate calcu- 
lating expression by using as weights said second 
membership grades as determined; and 

means for deriving a weighted mean of said weighted 
flow rate mean values calculated in accordance 
with said mean flow rate calculating expression for 
said flow velocity levels by using as a weighted 
said first membership grade as determined. 

23. A flow rate measuring system according to claim 
21, wherein said plurality of detecting means are dis- 
posed in a same plane extending orthogonally to the 
longitudinal axis of said pipe. 

24. A flow rate measuring system according to claim 
21, wherein said plurality of detecting means are dis- 
posed along a direction in which said fluid flows within 
said pipe. 
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